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Abstract. 

In this paper we consider the problem of a steady MHD flow of a non-Newtonian power-law and 
electrically conducting fluid in presence of an applied magnetic field. The boundary layer equa- 
tions are solved in similarity form via the Lyapunov energy method, we show that this problem has 
an infinite number of positive global solutions. 
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1. Introduction 

The study of non-Newtonian fluid flows has considerable interests, this is primarily because of the 
numerous applications in several engineering fields. Such processes are wire drawing, glass fiber 
and paper production, crystal growing, drawing of plastic sheets etc. For more details about the 
behavior in both steady and unsteady flow situations, together with mathematical models, we refer 
the reader to the books [[TJ by Astarita and Marucci, J21 by Bohme and the references therein. One 
particular non-Newtonian model which has been widely studied is the Ostwald-de Wael power-law 
model SHI, which relies the shear stress to the strain rate u y by the expression 

T~xy = k\lly\ Uy, (1-1) 

where A; is a positive constant, and n > is called the power-law index . The case n < 1 is referred 
to pseudo-plastic or shear-thinning fluid, the case n > 1 is known as dilatant or shear-thickening 
fluid. The Newtonian fluid is a special case where the power-law index n is equal to one. In the 
present work we shall restrict our study to the case n > 1. 

The magnetohydrodynamics (MHD) flow problems find also applications in a large variety of 
physical, geophysical and industrial fields [5]. It is also interesting to study the flow of non- 
Newtonian fluids with externally imposed magnetic fields. To the author knowledge MHD flow 
of non-Newtonian fluids was first studied by Sarpkaya 0. In Sapunkov derived the equations 
describing the similarity solutions for the non-Newtonian flow when the external applied magnetic 
field varies as x 12 ^ - , in presence of a pressure gradient, he used the method of series expansion. 
Later, Djuvic [8 J employed a Crocco's variables to study the unsteady flow with exponentially ex- 
ternal velocity (in time). Recently, Liao [9| introduced a powerful technique (homotopy analysis) 
to give analytic solutions of MHD viscous flows of non-Newtonian fluids over a stretching sheet. 
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In this paper, we reconsider the steady two-dimensional laminar flow of an incompressible viscous 
electrically conducting dilatant fluid over a stretching flat plate with a power-law velocity distribu- 
tion in the presence of a perpendicular magnetic field. Our interest in this work has been motivated 
by the work of Chiam ifTOl . who have considered the flow over an impermeable flat plate, for which 
similarity solutions were found via the Crocco transformation. 

2. Derivation of the model 

Consider a steady two-dimensional laminar flow of an incompressible dilatant and electrically con- 
ducting fluid of density p, past a semi-infinite flat plate. Let (x, y) be the Cartesian coordinates of 
any point in the flow domain, where x— axis is along the plate and y— axis is normal to it. Assume 
that a magnetic field H(x), is applied normally to the plate. 

The continuity and momentum equations can be simplified, within the boundary-layer approxima- 
tion, into the following equations (see fTllfTOl ) 

u x + v y = 0, (2.1) 

aji 2 H 2 

uu x + vu y = u(\u y \ n ~ u y )y + u e u ex H (u e — u) . (2.2) 

Accompanied by the boundary conditions 

u(x, 0) = U w (x), v(x,0) = V w (x) and u(x, y) — > u e (x) as y — > oo. (2.3) 

Where the functions u and v are the velocity components in the x and y directions respectively, 
u e (x) = UqcX™ is the free-stream velocity. The parameters u, n, fi, a and H are the kinematic 
viscosity, the flow behavior index, the magnetic permeability, the electrical conductivity of the 
fluid, and the magnetic field intensity respectively. The functions U w (x) = u w x m (u w > 0) and 

m(2n-l) — n 

V w (x) = v w x ™+i are the stretching and the suction/injection velocities respectively. 

In term of the stream-function which satisfied u(x, y) = ip y (x, y) and v(x, y) = —ijj x (x, y)), 

equations (12. II) , (12.21) can be reduced to the single equation 

Ipylpxy ~ Ipxlpyy = I '{\4 ^yyT' H >yy] ly + U e U ex H (u e - lf> y ), (2.4) 



P 



subject to 



m(2n — 1) —n 

i/j y (x,0) = u w x m , ifr x (x,0) = —v w x ~ L and ip y (x,y) = Uoox" 1 as y — ► oo. (2.5) 



According to Sapunkov 0, similarity solutions for problem (12.4I) . (|2.5I) exist only if the magnetic 

field has the following form H(x) ~ x^~ . 

To look for similarity solutions we define the following 

r}:=Ayx~ a and ijj(x, y) := Bx b f(r]), (2.6) 

where / is the transformed dimensionless stream function and rj is the similarity variable. Thanks 
to (12.61) . the function / satisfies the new boundary value problem 

(l/T-77 + aff" + m(l - f' 2 ) + M(l - f) = 0, 

(2.7) 

/(0) = a, /'(0) = 5, /'(oo) = 1, 



if and only if 



l + ra(2n-l) , l + m(n-2) 

a= , 6= , a — o = m, 

n + 1 n + 1 

and the parameters A and B satisfy 

AB = Uoo and uB^A 2 ^ = I. 

Where the primes denote differentiation with respect to rj, the function f'(rj) denotes the normal- 
ized velocity and the parameters 

afx 2 H 2 (n + 1) (n + l)v w u w 
M = , a = r- and 8 = — , 



u ooP (m + l)(pu 2 £- 



are respectively: The Hartmann number, the suction/injection and the stretching parameters. 
Such problems have been investigated by several authors for example, Anderson et al. ifTTTl . Zhang 
et al. [12] and Kumari and Nath lfT3l 

In the same context, Chiam [10J studied Problem (I2.1I) - (I2.3I) . To look for similarity solutions, he 
solved the following boundary value problem 

n \f»\n-i r + ff ' + p(i _ + M(1 _ //) = o, 

(2.8) 

/(0) = 0, /'(0) = 0, f(oo) = 0. 

Where /3 = (2 ^jff +1 . We aim here to stress that for n ^ 1, equation (I2.7l) i can be degenerate at 
some point rj s for which f"(rj s ) = (for more details see 03]]) and then any solution of (12.71) is not 
necessarily of C 3 (0, oo). Hence equations (I2.7l) i and (I2.8l) i are not equivalent. 
Let us notice that for the Newtonian case (n — 1), problem (12.71 ) reduces to the Falkner-Skan flow 
in Magnetohydrodynamics, which has been studied by Hildyard IfTTl , Aly et al. lfT8l and Hoernel 
|fP9l . The case m = M = leads to the generalized Blasius problem (see [20]). While the case 
m = —M, by a suitable scaling, is referred to the mixed convection of a non-Newtonian fluid in a 
porous medium (see for example [21 J). We note also that in absence of the magnetic field, problem 
(12.71) is simplified to the Falkner-Skan flow for non-Newtonian fluids. A complete study on this 
subject is given in [|22l by Denier and Dabrowski. 

Very recently, Aly et al. [18] reported a theoretical and numerical investigations on the existence 
of solutions to problem (12.71) for Newtonian fluids (n = 1), say 

/"' + ^ff" + m(l - f 2 ) + M(l - /') = 0, 

(2.9) 

/(0) = a>0, f'(0) = 5, f(0)= 7 - 
They showed that problem (12.91) has multiple solutions for any 5 E (0, V) and 76I satisfying 

2 777, 

~f 2 < —5 3 + M5 2 - 2(M + m)5, (2.10) 



where T = 



16/72 

1 h ^ 1 + 3M^ (m + M) 



> 1 . In the present work, we aim to extend their 



4m 

results to the non-Newtonian dilatant fluids (n > 1), by using a condition on 7 which is different 
from (12.101) and without any restriction on the parameter 6. 



3. Non-uniqueness of solutions 

Guided by the analysis of lfl4l . |[T5l and lfT6ll . we aim to prove the existence of solutions to problem 
(12.71) . for related values of the parameters m, M, n, a, 5 and 7. This result will be established by 
mean of the so-called shooting method, the boundary value problem (12.71) is then converted into 
the following initial value problem 

(l/T-77 + aff" + m(l - f' 2 ) + M(l - /') = 0, 

(3.1) 

/(0) = a, f'(0) = S, f'(0) = 7- 
Where the real number 7 is the shooting parameter. 

The initial value problem (13.11) can be transformed into the equivalent first order ordinary differen- 
tial system 

' /' = 9, 

< g >=\h\^rh, (3.2) 
k h' = -af\h\-m(l-g 2 )-M(l-g), 

with the conditions 

f(0) = a, g{0)=6, h{0) = \j\ n ~H (3.3) 

By the classical theory of ordinary differential equations, problem (13 .21) . (13 .31) has a unique local 
(maximal) solution for every 7 7^ 0. Let / 7 denotes this solution and (0, rj 7 ), r/ 7 < 00, denotes its 
maximal interval of existence. The main task now is to show how existence of solutions depends 
on 7. 

The local solution / 7 satisfies the following 

|r 7 r7" + a/y 7 - M(/ 7 + a) = |7r- x 7 + aa6 — (M + m) V + (a + m) P ' f^(r)dr. (3.4) 
Equation (|3.4I) will be used for proving the main results. 

Definition 3.1. A function / 7 is said to be a solution to OTTI) if f e C 2 (0, 00), \ f"\ n ~ l f" e 
C 1 (0, cx)) and satisfies 

lim /;(?7) = 1 (i) and lim /''(r/) = (ii) 
3.1. Suction/Injection flows (a G M) 

Theorem 3.1. Assume a e R, 5 > 0, M > , n > 1 and — ^ < m < -M. For any 7 
satisfying 

|7| n_1 7 > -a«5 (★), 
problem O.iD admits a global unbounded solution. 

Proof. From a physical point of view, it is more convenient to prove the result for the cases a > 
(suction) a < (injection) separately. 

We have to show that / 7 is a positive monotonic increasing function on (0, ry 7 ), globally defined 
and going to infinity with 77. For this sake we define the Lyapunov Energy function by 

1 m M 

V(v) = —-r\f"\ n+1 ~ ?/' 3 - + (M + m)f. (3.5) 



Which satisfies 

V\rj) = -aff" 2 . 

Then V is monotonic decreasing on (0, r^ 7 ). On the other hand, from equation (13.41 ) and condition 
(*) we see that / 7 ' and / 7 are positive on (0, 7/ 7 ) as long as / 7 exists. Using the Lyapunov function 
we see that _f 7 ' and / 7 " are bounded, since V is bounded from below by 3M + 4m . if J 7 were also 
bounded, say / L with L E (0, oo) (since / 7 is positive). Then / 7 '(r/) ^^oo which implies 
that f"(r)k) ~^koo 0, where (r] k ) k >o is a sequence tending to infinity with k. Using again (13.41) to 
deduce 

l/JfaOr 1 /^) + = -M(f y (vk) + a) + \i\ n - l i + aa6- 

(M + m)r] k + (a + m) ff(r)dr. 

Jo 

Letting k — > oo, the right hand side goes to zero while the left hand side goes to minus infinity, 
which is impossible, then / 7 is a global unbounded solution to (13.11) . 

From the above /' and /'' are bounded and /' is monotonic increasing on (r)i, oo), for rji large 
enough. Then there exists / > such that lim^oo f'Jrj) = /, and there exits a sequence (Cfc)fc, 
tending to infinity with k such that lim^oo /'' (Cfc) = 0. Making recourse to the Lyapunov function 

V we get lim^oo f"(( k ) = 0. 

Assume now that f" is not monotonic on any interval [772, 00). Then, there exists a sequence (r k )k 
going to infinity with k such that: 

• (l/?l^- 1 /?) / (^) = o, 

• l/ 7 T 1 f"( T 2k) is a local minimum, 

• |/ 7 T _1 / 7 '( r 2fe+i) is a local maximum. 
From (|3.1l) i. we have 

f/>( , m(l-/ 7 2 (r fc )) + M(l-/;(r fc )) 

Since / 7 is bounded and / 7 goes to infinity with 77^, we get easily from the above that f" goes to 
zero with 77. 

Now we show that / 7 satisfies (i). Recall that /' is a positive bounded function then /' — > I with 
I E (a, 00). At infinity we have / 7 ~ 77/ and from identity (13.41 ) we get 

l/ 7 T7" ~ ??H 2 + Ml - (M + m)) + o(l) 

as 77 approaches infinity, this leaves only the possibility that I is either 1 or — — — 1, thanks to the 
positivity of /' we deduce that 1 = 1. 

To finish we show the result for a < 0. In such case, the function f y is negative on a small 
neighborhood of zero. According to (|3.4I) / 7 cannot have a local maximum, then two possibilities 
arise: 

• Either / 7 < V77 G (0, ?7 7 ) 

• Or 3 77* such that / 7 < on (0,77^), f 7 (r]*) = and / 7 > V77 > 77^. 
Assume that the first assertion holds, then a > / 7 (oo) < and /'(oo) = 0, /' being positive we 
use again (13.41) to get that /'' is positive. A contradiction. Then, / 7 has exactly one zero 77^. We 
define the shifted function h by : 

77 1 — ► /i(?7) = / 7 (t7 + 77^), 



which satisfies /i(0) = 0, h'(0) = 5 and h"(0) > 0, and we use the above analysis to conclude that 
h is an unbounded slobal solution to (13.11). ■ 



3.2. Reversed flows (5 < 0) 



Now we pay attention to the case of reversed flows (5 < 0) . First, we show that the shooting 
parameter has to be positive. 

Proposition 3.1. Let / 7 be a solution to (fill) witfj m G (-^, -M), a < 0, 5 < and 7 < 0, 
dien the condition (i) is failed. 



Proof. Let 6 < 0, if 7 < then / 7 " is negative on some (0, 770), for r/ small, and equation (|3.1I) 
can be written as 

e F 

— \f* 
n 



(./■ 



,V7 



f//|l-n 
'7 I 



[m(i - /; 2 ) + m(i - /; 



a /' ?? 

where ^(77) = — f y \f y "\ dr. From this we see that 77 1 — > f"e F decreases and then f"(rf) < 

for all 77 G (0, 77 7 ). It follows that /' is decreasing on (0, ry 7 ) and then the condition (i) could not 
be satisfied. ■ 



Theorem 3.2. Let 5 < 0, a > and m G 

1 



3n' 



— M). For any j > satisfying 

M 



a-f - -5 2 7 n " 1 + aa 2 5 a 2 > (**), 

problem &3.1\) has a global unbounded solution. 

Proof. Let / 7 be the local solution of (13.11) . define the auxiliary function 



g(v) = f^V'T' 1 

which satisfies 

G'(rj) = -(m+M)/ 7 + 



1 f/2|_f//|n-l 1 „t2fl M 2 



1/7 I ^7 I 



+ </; - — 



(3.6) 



2a + m + 



n — l)a 



7-7 J-yl r> Jl /-v 



2n 2n - 7 . 7 r (i-/; a ) + M(i-/;) 

(3.7) 

and G(0) > 0. 

Since /' < 0, the function / 7 is negative on a small neighborhood of zero. Assume that there exists 
rji E (0, 00) such that 



/ 7 (t/)>0, /;<0 Vt7G[0,t^) and ^(771) 



0. 



Hence G is a monotonic nondecreasing function on (0, 771) and then G(r)i) < 0. Then G(rj) < 
for all 77 G (0, 771), in particular G(0) < 0, which is a contradiction with (**). Therefore we have : 

• Either f y > and / 7 < V77 > 

• Or 3r] 2 > : / 7 > 0, / 7 < V77 G (0,772), ^(772) = and ^7(772) is a local maximum. 
Assume that the first assertion holds, then / 7 has a finit limit at infinity, say L G (0, 00) and there 
exists a sequence (xk)k > tending to infinity with k such that f'(xk) g° es to zero at infinity. If /' 

is monotonic (resp. non-monotonic on any interval (77, 00)) we get /' goes to zero at infinity and 
then f"(5k) goes to zero at infinity for a sequence (5k)k>o going to infinity with k (resp. f!!(Sk) = 
and f~($k) goes to zero at infinity). Because G(0) < G(5k), we obtain a contradiction by taking 
the limit as k goes to infinity. 

Now, we claim that the function / 7 cannot have a local maximum after 772. Actually, assume there 
exists 773 > 772 such that ^(773) is a local maximum. At this point the function G takes a negative 
value and satisfies G(rj 3 ) > G(0) a contradiction. Since / 7 is monotonic increasing after 772 we 
deduce as the above that is a global solution. 

Next, we argue as in the proof of Theorem. 13 . 1 I to show that / 7 is unbounded at infinity and satisfies 
(i) and (ii). ■ 



3.3. Flow with large initial velocity (5 1) 

In this subsection, we construct asymptotic solutions to problem (13.11) when the real 5 is very large. 
Adopting the method used in [23J by Aly et al., we assume that such solutions can be written under 
the following form 

fiv) = V + Cg(t), where t = £ s r?, £ = 5 - 1 and r,set 
Then problem (13.11) reads 

£(r+2»)(n~2)(|0//|n-y/y + ar] ^r g f> + agg n _ ( 2m + M)^ r+S ^ - m(f = 0, 

g(o) = ar r , g'(0) = e~ (r+s) , g'(oo) = 0. 



(3.8) 



Setting r = and s = ensures that the highest derivative remains present in the resulting 
problem. 

As £ goes to infinity, we deduce 

(\g"\ n - 1 g"y + agg" + mg' 2 = J 

(3.9) 

0(0) = 0, g'(0) = 1, g'(oo) = 0. 

Problem (13.91) describes the steady free convection flow of a non-Newtonian power-law fluid over 
a stretching flat plate embedded in a porous medium. In [15], it was shown that for m E (— 0) 
any local solution g, whith positive values of r (r = g"(0)), is global and satisfies the following 
asymptotic behaviour 

l + m(2n-l) 

Cjf(t) ~ t l + m(n-2) ; aS ^ oo_ 

Consequently, a solution / for positive 7 and large 5 (if it exists), may have the following large 
//-behaviour 

|~ 1 m(n-l) + 2 " 

f(jj) ~ ?7 1 + ((5 — 1) l+m(n-2) 7 yl+m(n-2) 



4. Concluding remarks 

Based on the similarity transformation approach, the boundary layer equations for flows of purely 
viscous non-Newtonian dilatant and electrically conducting fluids are investigated. Using a shoot- 
ing argument, it is shown that the relevant problem admits an infinite number of solutions ( ll24ll [|25l ll26ll 
and [27]), this is due to the arbitrariness of the shooting parameter 7. From a physical point of view, 
we underline that 7 = f"(0) originates from the local skin friction coefficient Cf x , and the local 

(u x m ) 2 ~ n x n 

Reynolds number Re x = — ; via the the formula 

In conclusion, we may expect that the solutions determined above are physically acceptable. How- 
ever, only experiments are able to prove their physical existence. 
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